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Abstract. The vanishing ideal / of a subspace arrangement V\ U V2 U • • • U V m C V 
is an intersection Ii n I2 fl • • • fl I TO of linear ideals. We give a formula for the Hilbert 
polynomial of / if the subspaces meet transversally. We also give a formula for the 
Hilbert series of the product ideal J = ■ • ■ I m without any assumptions on the 
subspace arrangement. It turns out that the Hilbert series of J is a combinatorial 
invariant of the subspace arrangement: it only depends on the intersection lattice 
and the dimension function. The graded Betti numbers of J are determined by the 
Hilbert series, so they are combinatorial invariants as well. We will also apply our 
results to Generalized Principal Component Analysis (GPCA), a tool that is useful 
for computer vision and image processing. 

1. Introduction 

Suppose that V is an n-dimensional A-vector space. A subspace arrangement is a 
union 

A = V x u • • • u V m 

where Vi is a subspace of V for all i. Interestingly, various algebraic and topological 
properties of the arrangement A only depend on the dimensions rig := dim^ f) iGS Vi, S C 
{1,2, . . . , m}. Such properties are called combinatorial invariants of the subspace ar- 
rangement. For example, if K — R, then the topological Betti numbers of the comple- 
ment V\A are combinatorial invariants (see [Zj). If K = C, then the cohomology ring 
of V \ A is a combinatorial invariant (see fUJ |2]). For more on subspace arrangements 
and hyperplane arrangements, see [12 . 

Let Ij C K [V] be the vanishing ideal ^cy for j — 1, 2, . . . , m. The vanishing ideal 
/ of A is equal to the intersection Ji n h H • • • H I m . We also define J = 1% I2 • ■ ■ I m as 
the product ideal. We give a formula for the Hilbert series of J ( Theorem 14. Ij) . We also 
will give a formula for the Hilbert polynomial of I if all subspaces meet transversally 
( Theorem 14 .6J) . The Hilbert series of J is a combinatorial invariant ( Corollary 13. 2|) . 
but the Hilbert polynomial of / is not ( Example 15 .4|) . The Betti numbers and graded 
Betti numbers of J are also combinatorial invariants (Corollary I4.4|) . 

The original motivation for this paper comes from computer vision. A generalization 
of Principal Component Analysis naturally leads to the question of recovering the 
dimensions ni := dim Vi, i = 1, 2, . . . , m, given the Hilbert polynomial of the subspace 
arrangement. For more on Generalized Principal Component Analysis, see [To] . 
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2. Hilbert functions, series and polynomials 

Suppose that V is an n-dimensional vector space over a field K. We identify the 
coordinate ring R := K\V) with the polynomial ring := K[x\, X2, ■ ■ ■ , x n ] in n 
variables by choosing a basis in V. There is a natural grading R = © deN Rd where N = 
{0, 1, 2, . . . } is the set of natural numbers and Rd denotes the space of homogeneous 
polynomials of degree d. Let Z be the integers and suppose that M = © deZ Md is a 
finitely generated graded i?-module. We have Md = for d <C because M is finitely 
generated. The Hilbert function Hm of M is 

h M (d) = h(M, d) = dim K M d , d G Z. 

The Hilbert series of M is defined by 

H(M,t) :=^h(M,d)t d . 

d&L 

It is a Laurent series because h(M, d) = for d C 0. 

Example 2.1. For M = R we get h(R,d) = ( n ~l\ d ) for d > and h(R,d) = for 
d < 0. So we have 

00 1 

«w*) = E("^«'=7rhvr- 

Define the polynomial Iir by 

~ ~ (rf + 7l-l)(rf + 7l-2)---(d+l) 

ftji(a) = h(R,d) = — . 

[n — 1)1 

One easily checks that hn(d) = hft(d) for d > 1 — n. 

Let M be again a finitely generated graded i?-module. For r £ Z we define the 
shifted module M[r] = dgZ M[r] d by M[r} d := M r+(J , d G Z. Shifting the degrees 
affects the Hilbert function and the Hilbert series as follows: 

h(M[r],d) = h{M,d + r), d,reZ, 

H(M[r],t) =r r H(M,t). 
The module M has a minimal finite free graded resolution 

(1) 0^0 R[-j}^ - - • • • - R[-j]*u - M - 0. 

by Hilbert's Syzygy Theorem (see for example [E],fEl §13], |3 §19.2]). The nonnegative 
integers (3^ are called the graded Betti numbers. For all but finitely many pairs 
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we have flij = 0. The Betti numbers are defined by $ = X^ezAj ( n °t to be confused 
with the topological Betti numbers of the complement of the subspace arrangement 
mentioned earlier). Without loss of generality we may assume that f3 r 7^ 0. The 
nonnegative integer cd(M) := r is the cohomological dimension of the module M and 
is at most n. The Castelnuovo-Mumford regularity of M is 

reg(M) := max{j — i | < i < r, fcj ^ 0}. 

From the exactness of (0) follows that 

r r 

(2) h(M, d) = ^(-i) 4 PuKR[-M = J2^y E rf - i) 

i=o jez i=o jgz 

and 



TC(M,t) 



{i-ty 

We define the Hilbert polynomial Hm of M by 

r 

(3) h M (d) = h(M, d) = ^(-l) 1 d-j). 

i=o jez 

Corollary 2.2. If d > 1 — n + reg(M) + cd(M), then we have 

h M (d) = h M (d). 

Proof. If fcj ^ 0, then 

d-j = d-i-(j-i)>d- cd(M) - reg(M) > 1 - n, 
hence h^(d — j) = hn(d — j) (see Example 15 .lp . The corollary follows from @ and 

3. Subspace arrangements 

For the remainder of this paper, let V be an n-dimensional vector space and suppose 
that Vi, V 2 , . . . , V m are subspaces of V. For a subset X CV, let T(X) C R = K[V] be 
its vanishing ideal. Define Ij = 1{Vj) for j = 1, 2, . . . , m. The union 

A = V x U V 2 U • • • U V m 

is a subspace arrangement. Its vanishing ideal is 

i:=x(A) = Jin J 2 n •••ni m . 

Define 

J := ■ • • Jm- 

Theorem 3.1. TVie regularity of I and J are bounded by reg(J) < m cm<i reg(J) < m. 
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For reg(J) < m, see 11} . The bound reg(J) < m was proven in jBlUj- 

For a ideal a C i? we have cd(o) = cd(i2/a) — 1 < n — 1. In particular, we get 
cd(/) < n — 1 and cd(J) < n — 1. 

Corollary 3.2. We /mve = /i/(c/) and hj(d) = hj(d) for d > m. 

Proof. By Corollary I2~21 hj(d) = hj(d) and hj(d) = hj(d) for 

d > m = (1 — + m + — 1) > 1 — n + reg(I) + cd(J). 



□ 



4. Main results 



For 5 C {1,2, ... ,m}, define Is = f] seS I s and J s = ELes^- Note that / = 
I{i,2,...,m} an d J = J{i,2,...,m}- We use the convention Iq, = J® = R. For S C {1,2,..., m} 
define = Hies — dim Vs and eg = n — n s is the codimension of in V. We 

also set n, = n{i} = dim V* and q = C{i} = n — n« for i = 1, 2, . . . , m. 

We define polynomials ps(t) recursively as follows. First we define 

m{t) = i- 

If S 1 7^ and px(^) is already defined for all proper subsets X C S, then ps(£) is 
uniquely determined by 

£(-t)'*W(*) = mod (1 - t) Cs , deg(p s (t)) < c s . 

xcs 1 

Here deg(px(^)) is the degree of the polynomial pxif) and |X| is the cardinality of the 
set X. 

Theorem 4.1. We have 



where p(t) =P{i, 2l ...,m}(*)- 

Corollary 4.2. T7te Hilbert series TC(J,t) depends only on the numbers n$,S C 
{1,2,. ..,m}. 

Proof. This follows immediately from the recursive formulas for ps{t), S C {1, 2, . . . , n} 
and the observations that n = n% and cs = n — n$ for all 5 1 . □ 

Proposition 4.3. Let {Aj} fre £ne graded Betti numbers of the ideal J. If faj ^ 
£/ien j = m + i and fa := . faj = fa, m +i is the i-th Betti number. We have 

„, r A _ EI =0 (-l)^ +m _ t m p(t) 
n{,) (l-t)» (l-t)"' 

where 

P {t) =po- fat + p 2 t 2 ---- + {-i) r faf, 



HILBERT SERIES OF SUBSPACE ARRANGEMENTS 5 

andr = cd(J) < n — 1. The Betti numbers and graded Betti numbers of J are uniquely 
determined by the Hilbert series of J, so they are combinatorial invariants. 

Proof. By Theorem 13.11 we have 

reg(J) < m. 

The ideal J = Ji J2 ■ ■ ■ I m is generated in degree m. This means that J has a linear 
minimal free resolution (see [SJ Proposition]): 

R[-m - rf r -> ► R[-m - if 1 -> i^-mp -> J -> 0. 

The proposition follows. □ 

Corollary 4.4. JTie J?e££z numbers and graded Betti numbers of J are uniquely deter- 
mined by the Hilbert series of J, so they are combinatorial invariants. 

Definition 4.5. The subspaces V\, . . . , V m are called transversal if 

c s = min (n, J2 ieS 
for all S 1 C {1, 2, . . . , m}, where min denotes the minimum. 

Note that we always have cs < min (n, ^ ig5 cA ■ So the subspaces are transversal if 
any intersection of some of the subspaces has the smallest possible dimension. 

Theorem 4.6. Suppose that Vi,...,V m are transversal. Then Tt(I,t) — fit) and 
TC(J,t) — f(t) are polynomials in t, where 

f(f) _ nti(i -a-*)*) 

Corollary 4.7. IfV\,..., V m are transversal, then 

hj(d) = hj(d) = hj(d) = hj(d) 

for all d > m. 

Proof. From Theorem 14.61 and the exact sequence 

0^ J -> / -> Ij J ^0 

follows that 7i(I / J, t) = 7~C(I, t)—7i(J, t) is a polynomial. So J/J has a finite dimension. 
It follows that hi(t) = hj(t). We have 

hi(d) = hj(d) = hj(d) = hj(d) 
for d > (1 — n) + m + (n — 1) = m by Corollary 13.21 □ 

Corollary 4.8. IfVi,..., V m are transversal, then 

where cs = X^es Ci an< ^ ^ e sum ^ s over a ^ ^ C {1, 2, . . . , m} for which cs < n. 
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5. Examples 

Example 5.1. Suppose that m = 3 and n = 3 and that Vi, V2, V3 are one-dimensional 
and distinct. Now Vi, V2, V3 correspond to 3 points Pi,P 2 ,P 3 G P 2 . Suppose that 
Pi, P2, P 3 are not colinear. We can change coordinates such that 

P 1 = (1 : : 0), P 2 = (0 : 1 : 0), P 3 = (0 : : 1). 

We have 

7 = ii n h n i 3 = (#2, £3) n (xi, a? 3 ) n (xi, x 2 ) = (x 2 x 3 , Xi2 3 , 2ix 2 ) C if[xi, x 2 , x 3 ] 
and 

J = IlI 2 h = (xiX2X3,xjx2,xlx3,xlx3,xlx 1 ,xlx 1 ,xlx 2 ) C. K[xi,X 2 , X S ]. 

The minimal free resolutions follows 

-> P[-3] 2 -> P[-2] 3 -> / -> 0, 

-> P[-5] 3 -> P[-4] 9 -> P[-3] 7 -> J -> 0. 
So we have cd(I) = 1, cd(J) = 2, reg(i) = 2, reg(J) = 3, 

, r N 3i 2 -2t 3 T . 7t 3 -9i 4 + 3t 5 

w(/ ' f) ^ 0^ andW(J - f)= (i-t)' • 

The Hilbert functions, which can be found as the coefficients of the Hilbert series, are 
given in the following table. 



d 





1 


2 


3 


4 


5 


hi(d) 








3 


7 


12 


18 


hj(d) 











7 


12 


18 



The Hilbert polynomials hj(d) and hj(d) are both equal to 

d 2 + 3d - 4 



Let us verify Theorems 14.11 and 14.61 The subspaces V\ , V2 , V 3 are transversal. Let 



w /(*)=nSi 



(l-(l-t) c (l-(l-t) 2 ) 3 -2 + 6t-3t 



-t)> 



(i-ty 



+ 2 - 3t 2 + t 3 . 



We can write 



and 



„„ T , 3t 2 - 2t 3 -2 + 6t-3t 2 
W(J, t) = — = - ^— + 2 



nun 7/: ': {)/ '.: :</r ' 2 .; (i/ .:;^-2 



(l-t)3 

As Theorem 14.61 predicts. 7i(I,t) — fit) and TC(J,t) — fit) are polynomials. 
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Using the recursive definitions of ps we get p$(t) = 1, 
P$(t)~tp{i}(t) = l-tp {2} (t) = l-tp {3} (t) = mod (1-t) 2 , deg(p { i } ),deg(p {2 }),deg(p {3} ) < 2, 
so P{iy(t) = P{2}(t) = P{3}(t) = 2 — t. Furthermore, 

1 - t(p {1} (t) + p {2} (t)) + t 2 p {lt2] {t) = mod (1 - t) 3 , 

so we have 

t 2 p {ly2 }(t) = -1 + 2t(2 - t) = -1 + 4t - 2t 2 mod (1 - t) 3 
and deg(p{i j2 }) < 3. It follows that 

P{1,2}(*) = P{1,3}(*) = P{2,3}(*) = 4 - 4t + t 2 

Finally we have 

l-t(p {1} (t)+p {2} (t)+ m (t)H^^ = mod (1-t) 3 , 

so 

^ 3 P{i,2,3}(t) = 1 - 3t(2 - t) + 3t 2 (4 - 4* + t 2 ) = -2 + 6* - 3t 2 mod (1 - t) 3 . 
From this follows that 

P{i,2,3}0) = 7-9t + 3t 2 . 

Theorem 14.11 correctly gives 

„,j +\ _ * 3 P{i,2, 3 }(*) _ 7t 3 - 9t 4 + 3t 5 

"■^V- (l_ t )3 - (l_ t )3 • 

Example 5.2. If Pi,P 2 ,P 3 are colinear, then, after a change of coordinates, we may 
assume 

P x = (1:0,0), P 2 = (0:1:0), P 3 = (1 : 1 : 0). 

We have 

I = Iin/ 2 n/ 2 = (x 2 ,x 3 ) n (a;i,x 3 ) n (x x -a; 2 ,x 3 ) = {x 1 x 2 {x 1 -x 2 ), ar 3 ) C A"[xi, x 2 , x 3 ] 
and 

J = JiJ 2 / 3 = (xix 2 (xi - X2),x 2 x 3 ,XiX 2 x 3 ,X2a: 3 ,a:iX 2 ,a: 2 a: 2 ,X3) C x 2 , x 3 ]. 

We have minimal free resolutions 

-> P[-4] -> 1] © P[-3] -»• I -»• 0, 

-> P[-5] 3 -> P[-4] 6 -> P[-3] 7 -> J -> 0. 
We get cd(I) = 1, cd(J) = 2, reg(J) = reg(J) = 3, 

. . t + t 3 -t 4 ,„,. T . 7t 3 -9t 4 + 3t 5 
H(I, t) = (1 _ t)3 and W( J, t) = (1 _ t)3 • 
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The Hilbert functions are 



d 





1 


2 


3 


4 


5 


hj(d) 





1 


3 


7 


12 


18 


hj{d) 











7 


12 


18 



We can compute fit) and it is the same as in Example 15.11 ®- Now 

t + t 3 - t 4 -2 + 6t - 3t 2 



H(I,t) 



+ 2 + t, 



(i-t)* (i-ty 

so 7Y(7, t) — f(t) is indeed a polynomial as in Theorem 14.61 Theorem 14. II gives the same 
result for TC(J,t) as in Example [T] 

Example 5.3. Let n = 4 and m = 3. Consider the ideals I\ = (x2,X3), I2 = (xi,x^) 
and ^3 = (xi,X2) in the polynomial ring K[x\, x%, x%, X4] in 4 variables and let again 
/ = l\ fl li fl J 3 . The ideals Ji, J 2 , -^3 correspond to 2-dimensional subspaces Vi, V2, V 3 of 
the 4-dimensional space V such that Vifl V 2 n V 3 is 1-dimensional, and V\ + V2 + V3 = V. 
We get similar free resolutions of / as in Example 15.11 

-> i?[-3] 2 -> i?[-2] 3 -»• I -»• 0, 



We have 



H(I,t) 



3t 2 - 2t 3 



Example 5.4. Let n = 4 and m — 3. Consider the ideals /1 = (xi,x^), I2 = (£2,^3) 
and ^3 = {x\ — X2,xz) in the polynomial ring K[x\,X2,x^x^ in 4 variables and let 
I — Ii fl I2 D J3, J = hhl-i- The ideals ii, J 2 , J 3 correspond to 2-dimensional subspaces 
^1)^2,^3 of the 4-dimensional space V such that Vi fl V2 fl V3 is 1-dimensional, and 
V\ + V 2 + Vs is 3-dimensional. In projective space P 3 , we have 3 lines lying an a plane 
and going through 1 point. We get similar free resolutions of / and J as in Example 15 .21 

_> R[- A ] _> R[-l] 

R[-hf -> i?[-4] 9 











-3]- 

-3] 7 



-0, 
J -> 0. 



We obtain 



W(J,t) = 
Taking the difference gives 



£ 4 



(l-f)< 



and W(J,t) 



7t 3 - 9t 4 + 3i 5 



(l-f)< 



H(I,t) -H(J,t) 



3t 2 + t 



which is not a polynomial. Note that Theorem 14. (il does not apply because Vi, V2, V3 
are not transversal. We have 

«('.«> = %5« 



t 4 
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but in Example 15. HI we got 

«('.«) = 

This shows that Ti{I, t) is not a combinatorial invariant. The difference of both Hilbert 
series is not even a polynomial. This implies that the Hilbert polynomial hi(d) is not 
a combinatorial invariant either. 

6. Complexes of product ideals and intersection ideals 
Theorem 6.1 (See Chapter IV of ^3]). There exists complexes 

°^ J ^ Is - • • • - Is - R - 

|S|=m-l |S|=m-2 |S|=1 

and 

|5|=m-l |S|=m-2 |S|=1 

whose homologies are killed by a = YlJLi Ij- 

To describe the the maps in the complexes in Theorem 16 .11 it suffices to define maps 
It — ► Is and J T — > J5 for all subsets 5, T C {1,2, ...,m} with |T| = \S\ + 1. If 
^ = {ii, « 2 , • • • , i r } with zi < i 2 < ■ ■ ■ < i r and S = i 2 , . ■ . , i s -i, i s +i, ■ ■ ■ , i r } 
then the maps It — > -Ts and — > in the complexes in Theorem 16.11 are given by 
/ i— > (— l) s /. All other maps are equal to 0. 

Corollary 6.2. If\) : -, m] = f)7=i V i = ^en 

£ (-lf\H(I s ,t) 

SC{l,2,...,m} 



and 

are polynomials in t. 
Proof. The ideal 



^ (-l)'^(J 5 ,t) 

SC{l,2,...,m} 



E^= J (nr=i^)=^(W) = m 

i=l 

is the maximal homogeneous ideal of R. 
Suppose that 

is a complex of finitely generated graded .R-modules. The i-th homology group is 

Hi = ker(<9;)/im(<9 m ). 
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We have 

r r 

2=0 1=0 

If raHi = 0, then Hi is finite dimensional, and Ti,(Hi,t) is a polynomial for all i. 
We now apply this to the complexes in Theorem 16. II □ 

Proof of Theorem \4.1\ By Proposition 14.31 we can write 

EI=o(-i)W +m _ *>(*) 



H(J,t) 



(l-t) n (l-t)™' 
where 

(5) P (t) =p Q -p 1 t + --- + (-l) r f3rf 

is a polynomial of degree r < cd(J) < n — 1. Similarly we can write 
with 

(6) deg(p s (t)) < n - 1 

for all 3 C {1,2,..., m}. 

Let W = V/V s and define W x = V x /V s for_X C 5. Let J { C = Hft"-^] be 

the vanishing ideal of for all i E S. Define J x = ELex ^ ^ or a ^ X C. S. 

We have 

inside = ® From this follows that 

t W Px(t) _ ajfT ^_ H(J x ,t) 

In particular, we have 

n(J s ,t) - (1 _ t)n _ ns ■ 

From this it follows that deg(ps(t)) < dimW — 1 = n — n$ — 1 (see (jHJ)). Since 
P|j G5 = in W, Corollary IO implies that 

xcs xcs \ ' 

is a polynomial in t. Multiplying with (1 — t) n ~ ns gives 

^2i-t) W Px(t) = mod (1 -t) n ~ ns . 

xcs 

□ 
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Proof of Theorem \4.6\ 

Special case: Suppose that c\ + • ■ ■ + c m < n. After a change of coordinates (and 
using that the arrangement is transversal) we may identify with 

%[*] ,g K [C2] ® • • • ® K [Cm] ® K [n - Cl —- Cm] 

and Ik with 

if [ci1 <g> • • • ® K [Ck - l] <g> m Cfc «) iT [Cfc+l] <g> • • • ® if [Cm] ® if 

for /c = 1, 2, . . . , m. Here m r is the homogeneous maximal ideal of K M. We get 

J = m Cl (8) m C2 <g> • • • ® m Cm ® K [n ~ Cl Cm] = J. 

Note that H(K^,t) = (1 - fT r and H(m r ,t) = (1 - fp r - 1 for all r. Therefore, we 
get 

W(J,t) =H(J,t) = H(m Cl ,t)H(m C2 ,t) ■ ■ ■H(m Cm ,t)H(K^~---- c -\t) = 

n ((! - - 1) V a - t) ci+ - + — - = Ei( ( 1 1 : ( ^~ t)c ' ) - 

The general case. We prove the Theorem by induction on m. The base case m — 1 
follows from the special case above. If c\ + • — h c m < n then we are also done. Let us 
assume that c\ + ■ • • + c m > n. In particular we have C[i t 2,...,m} — min(n, Y^T=i c «) = n i 
so flUi = (0)- By Corollary 16.21 we have that 

(-lf l H(Is,t) 

SC{l,2,...,m} 

is a polynomial. 

By induction we have that 

Hi s (t) -wJJ^Vi 

is a polynomial for all strict subsets S C {1,2,..., m}, where Uj = 1 — (1 — t) c,; and 
w = (1 — t) _n . To show that 

m 
t=l 

is a polynomial, it suffices to show that 

(7) e (-i) isi w/ s ,t)-t«n^) 

SC{l,2,...,m} «GS 

is a polynomial. 

Now 

E (-i) |5| ^(^,t) 

SC{l,2,..,m} 
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is a polynomial by Corollary E21 and 

m m 

j2 (-i) 151 ^ = w n^ 1 - v *) f y n n^ 1 - *)* = (* - tr + - +c -- n 

SC{l,2,...,m} i&S t=l i=l 

is a polynomial as well. Therefore (J7J) is a polynomial. □ 

7. Application to Generalized Pricipal Component Analysis 

The object of Principal Component Analysis (PCA) is to approximate a data set 
inside a vector space V by a subspace of smaller dimension. In Generalized Princi- 
pal Component Analysis (GPCA) one tries to approximate a data set inside a vector 
space V by a union of subspaces spaces (in other words, a subspace arrangement). 
Some applications of GPCA are motion segmentation (see [T7| image segmenta- 
tion (see |14j). image compression (see 9J) and hybrid control systems For an 
overview of GPCA, see [To] . 

A first start in GPCA is to decide on the number of subspaces and the dimensions 
of the subspaces of the subspaces arrangement that will approximate the data. 

Suppose that V\, v%, . . . , v r G V are data points. Here r is fairly large. Suppose that 
Vi, . . . , v r are contained in some subspace arrangement A — V\ U • • • U V m , unknown to 
us. We would like to recover m, . . . ,n m where rij = dimV^. Let <Xj be the vanishing 
ideal of the ray through Vj. Then we have that 

h(a! n • • • n a r ,d) = h(I,d) 

for small values of d, where / = 1{A) as before. Now 

h(a.\ n • • • n o r , d) 

can be computed using linear algebra for small values of d. Therefore, we can determine 
h(I, d) for small values of d. So an important question is, given h(I, d) for small 
values of d, can we determine the dimensions ni,n 2 , ■ ■ ■ ,n m 7 Proposition 17. II gives an 
affirmative answer if the subspaces are transversal. Of course, in real applications the 
data is approximated by the subspaces arrangement, but not contained in it. In that 
case, using the PCA method in K[V]d one can still can estimate the value h(T(A), d). 

Proposition 7.1. Assume that the arrangement is transversal. Suppose that ci, . . . , c m 
are unknown, but we know the values of the Hilbert polynomial 

h^d) 

for d — m, m + 1, . . . , m + n — 1, then we can recover ci, . . . , c m . 

Proof. Note that hj(d) = hj(d) for d > m. Since we know hj(d) for d = m,m + 
1, . . . , m + n — 1 and hi has degree < n — 1, hi is uniquely determined. From this, we 
can determine TC(I, t), up to a polynomial. Suppose that t) is equal to a(t)/ (1— t) n 
up to a polynomial. Let b(t) be the remainder of division of a(t) by (1 — t) n . Then b(t) 
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has degree < n and Ti(I, t) is equal to b(t)/ (1 — t) n modulo a polynomial. So we have 
that 

d 

w) = n^ 1 - ( x - *) ci ) m ° d ( x - *) n 
i=i 

and 

6(1 - t) = - t Cl ) mod t n . 

1=1 

Let rj be the number of the c/s equal to i Then we have 

d 

6(1 - i) ee - t ) ri mod t n . 

i=i 

From this we can easily determine ri, r 2 , . . . , r n _i (in that order). Indeed, if we already 
know ri, . . . , r s , then the Taylor series of 

&(1 ~ t) 

is 

1 — r s+ it s+1 + higher order terms. 
So we find the value of r s+1 . □ 
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